First Results on Double Beta Decay Modes of Cd, Te and Zn Isotopes with the 

COBRA Experiment 



O 

o 



00 



X 



T. Bloxham,! A. Boston,"^ J. Dawson,^ D. Dobos,^ S.P. Fox/ M. Freer, ^ B.R. Fulton/ C. Gofiling/ 
P.F. Harrison, 6 M. Junker, ^ H. Kiel,^ J. McGrath,'^ B. Morgan,'^ D. Munsterniann,^ P. Nolan," 
S. Oehl,2 Y. Ramachers,^ C. Reeve,'^ D. Stewart,^ R. Wadsworth,'^ J.R. Wilson,^ and K. Zuber^ 

^School of Physics and Astronomy, University of Birmingham, B15 2TT, UK 
^ Lehrstuhl fiir Experimentelle Physik IV, Universitdt Dortmund, Otto-Hahn Str. 4,44^27 Dortmund, Germany 
'Laboratori Nazionali del Gran Sasso, S.S. 17 BIS km. 18.910, 67010, Assergi, L'Aquila, Italy 
^ Dept. of Physics, University of Liverpool, Liverpool L69 7ZE, UK 
^ Dept. of Physics and Astronomy, University of Sussex, Brighton, BNl 9QII, UK 
^Dept. of Physics, University of Warwick, Coventry CV4 7AL, UK 
''Dept. of Physics, University of York, Heslington, York, YOlO 5DD, UK 
(Dated: February 5, 2008) 

Four 1 cm'^ CdZuTe semiconductor detectors were operated in the Gran Sasso National Laboratory 
to explore the feasibility of such devices for double beta decay searches as proposed for the COBRA 
experiment. The research involved background studies accompanied by measurements of energy 
resolution performed at the surface. Energy resolutions sufficient to reduce the contribution of two- 
neutrino double beta decay events to a negligible level for a large scale experiment have already been 
achieved and further improvements are expected. Using activity measurements of contaminants in all 
construction materials a background model was developed with the help of Monte Carlo simulations 
and major background sources were identified. A total exposure of 4.34 kg-days of underground data 
has been accumulated allowing a search for neutrinoless double beta decay modes of seven isotopes 
found in CdZnTe. Half-life limits (90% C.L.) are presented for decays to ground and excited states. 
Four improved lower limits have been obtained, including zero neutrino double electron capture 



transitions of 
respectively. 

PACS numbers: 



Zn and "Te to the ground state, which are 1.19 x 10 years and 2.68 x 10 years 
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I. INTRODUCTION 



In recent years, a range of neutrino oscillation experi- 
ments [1, 2, j|] have successfully proved that neutrinos are 
massive particles. Although such experiments are sensi- 
tive to a mass-difference rather than absolute neutrino 
mass, the data suggest a neutrino mass eigenstate of at 
least 50meV. To further probe the neutrino's properties, 
it is necessary to look to other processes such as neutri- 
noless double beta decay (Ov^P), which violates lepton 
number by two units; observation of this process would 
confirm the Majorana nature of the neutrino and the rate 
of this rare decay is proportional to the absolute neutrino 
mass scale. For recent reviews of double beta decay see 
Refs. Hi and [|. 

The COBRA experiment uses CdZnTe (CZT) semicon- 
ductors to search for 0i//3/3 @ . CZT contains nine double 
beta emitters, five of which can decay via double beta de- 
cay, i.e. emitting two electrons, and four of them via ei- 
ther double electron capture, a combination of a positron 
emission with electron capture or double positron emis- 
sion. The study of the positron and electron capture 
modes can be used for lepton number violating decay 
searches on an equivalent level to 0^/3(3, however the 
phase space for the positron modes is strongly reduced, 
making them less sensitive. Nevertheless, is has been 
shown that the positron/electron capture mixed modes 
have an enhanced sensitivity to right-handed weak cur- 



rents and thus can help to disentangle the underlying 
physics mechanism of Oi^ (3(3 if observed Q- In addition, 
excited states transitions can be explored with high effi- 
ciency and low background using coincidence techniques 
among the detectors. These decays would allow an inde- 
pendent search for double beta decay searching for the 
de-excitation photons together with the electron signal. 
The main focus of the work described in this paper is the 
study of background through measurements performed 
underground, and energy resolution studies carried out 
in a surface laboratory. Optimisation of these quantities 
is vital for a successful search for Qv(3(3, because, in the 
background limited case, the observable half-life depends 
on them with a square root behaviour. In addition, half- 
life limits for seven double beta isotopes contained in nat- 
ural CZT have been determined from data collected with 
a small prototype detector accumulating an exposure of 
4.34 kg-days. 



II. SENSITIVITY 

Although COBRA is able to search for Ov(3(3 in a num- 
ber of isotopes, the sensitivity for the modes with lower 
Q-values will ultimately be limited by background con- 
tributions from two neutrino double beta {2v(3(3) decays 
of the isotopes with higher Q-values. The contribution of 
2v(3(3 decay events to the current data set is negligible, 
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FIG. 1: (Color online) The expected half-life sensitivity for 
Ot'/J/J-decay of '^^^Cd in an array of 64,000 1 cm^ detectors, 
enriched to 90% in ""^Cd for three scenarios of background 
and energy resolution {AE — FWIfM): 
A = 10~^ counts/(keV kg year), AE = 2% at 2.8 MeV. 
B = 10"^ counts/ (keV kg year), AE = 1% at 2.8 MeV. 
C = 5 X 10"" counts/(keV kg year), AE = 1% at 2.8 MeV. 
A half life of about 2x10^^ years corresponds to a neutrino 
mass sensitivity of about 50 meV using matrix elements from 
Ref. i). 



but for a very sensitive neutrino mass search, COBRA 
will focus on ^^^Cd, which has the highest Q- value of 
2809 keV for the nuclear decay to ""^Pd. A peak will 
occur at this energy in the sum energy spectrum for the 
case of Qv(3f3. 

Crucial experimental parameters, besides the mass of the 
detector/sample, are energy resolution and the number of 
contaminating background events in this range, as shown 
in Fig. [1] Possible background sources include cosmic 
rays, the natural radioactive decay chains, radioisotopes 
produced by cosmic ray interactions within the mate- 
rials used and neutrons. To get a first glimpse of the 
background using CZT detectors, a prototype setup has 
been installed in the Gran Sasso Underground Labora- 
tory (LNGS) in Italy, which provides an average shielding 
of ~3500mwe against cosmic ray sources. 



III. EXPERIMENTAL SETUP 

The data presented here were obtained with four 
Icmxlcmxlcm CZT semiconductor detectors, each of 
mass -6.5 g, provided by eV-PRODUCTS. They utilise 
coplanar grid technology to ensure that only the electron 
signal is read out ^] and hence symmetric energy peaks 
are obtained. All four crystals were operated at a voltage 
of — I500V with a 20-40 V grid bias applied between the 
two anodes. The crystal surfaces, except for the gold- 
coated cathode side, are covered with a passivation paint 
which prevents oxidation and deterioration in detector 
performance over time. 

The four detectors were mounted in a copper brick sep- 
arated from all preamplifier electronics by ~25 cm. The 



Subset Livetime (days) Events/(keV kg day) 
>500keV >600keV in 2-3 MeV range 



lA 


32.86 


150.99 


0.57 ± 0.02 


2A 





112.34 


0.59 ± 0.03 


3A 


52.84 


52.84 


0.38 ± 0.03 


4A 


62.62 


134.64 


0.34 ± 0.02 


IB 


16.03 


16.03 


0.54 ± 0.07 


IC 


197.53 


197.53 


0.56 ± 0.02 



TABLE L Livetimes for data sets prepared above two different 
energy thresholds (500 keV and 600 keV) , and an indication 
of the average background level in the 2-3 MeV region. 



copper brick was part of a (20 cm)"^ cube of electro- 
polished copper which was embedded in a further 15 cm 
of lead. The whole setup was located in a Faraday cage 
made from copper plates. The cage was surrounded by 
a neutron shield, consisting of 7 cm thick boron-loaded 
polyethylene plates and an additional 20 cm of paraffin 
wax at the bottom. This neutron shield was upgraded to 
cover 3 sides of the cage, as well as the base, with paraf- 
fin wax before the third data taking period (period C). 
Data collection commenced with a CAMAC based data 
acquisition (DAQ) system in which the signals were fed 
into four peak sensing ADC modules (LeCroy 3511 and 
3512) via custom-built preamplifiers and shaping main- 
amplifiers. Prior to period C, the system was upgraded 
to a VME-based DAQ with four custom-built, peak sens- 
ing, 14-bit ADC channels. 

IV. DATA ACQUISITION 

The data analysed in this paper can be divided into 
three periods, separated by upgrades to the experimen- 
tal configuration that could have affected the background 
contributions. In period A, the crystals were held in per- 
tinax mounting plates and connected directly to lemo ca- 
bles. These mounting materials were replaced by cleaner 
(in the radiopurity sense) delrin holders, whilst the lemo 
connections were exchanged for copper traces mounted 
on kapton foils at the start of period B. Before period C 
the data acquisition hardware was upgraded to the VME 
system and the paraffin neutron shielding was completed. 
During the latter periods not all crystals were fully op- 
erational, so only data from crystal 1 were analysed. In 
each period, individual runs were limited to one hour and 
selection criteria were applied on a run-by-run basis to 
reject data affected by "bursts" of abnormally high event 
rates. Subsequent studies have shown two main causes 
for such bursts: vibrations of the apparatus that cause a 
piezo-electric effect in the crystals resulting in false event 
signals, and breakdown effects due to faulty contacts to 
the crystal electrodes. To reject the affected data, firstly 
the dead-time was calculated from the total number of 
events per hour and the length of the event readout cycle 
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FIG. 2: (Color online) The prepared data sets for periods A, 
B and C. Period A comprises data from all 4 detectors and 
has slightly less events in the high energy regions as the back- 
ground levels in detectors 3 and 4 are lower. However, periods 
B and C show less events at low energies due to improvements 
in the cleanliness of mounting materials. 



during which new events could be missed. 207 runs in 
which the dead-time exceeded 2% (due to large numbers 
of events with low ADC counts) were rejected. For the 
remaining runs the distribution of the number of events 
per hour in the energy range 300-4000 keV was fitted 
with a Poissonian distribution. Runs with an event rate 
exceeding the 99% upper limit of the fitted distribution 
were discarded from the analysis. A further 203 runs 
were rejected in this way. No other cuts were applied to 
the data and in total 2.5% of the runs were rejected. 
It is possible that each crystal had different background 
contributions due to intrinsic and surface contaminants 
(as indicated by the values in table |T| . For this reason, 
data collected with each crystal, and in each period, were 
considered as a different data sub-set. The summed data 
for each of the collecting periods are shown in Fig. [51 

Although detectors have been operated for extensive 
periods with thresholds below 100 keV, at times it was 
necessary to raise the threshold for data collection to ex- 
clude electronic phenomena. Therefore, to maximise the 
livetime for analysis, two data sets have been prepared; 
one with an energy threshold of 600 keV, and one with 
a 500 keV threshold that omits any runs with thresholds 
in the range 500-600 keV. The livetimes for each subset 
and threshold are given in table [J The total livetime 
for the high-threshold data set is 4.34kg-days, whilst the 
low-threshold data set comprises 2.36kg-days. 
The energy resolution and stability of the detectors was 
calibrated regularly with the help of ^^^Cs, ^°Co and 
^^^Th sources. Time-averaged resolution functions were 
determined for each data subset with all crystals show- 
ing a linear increase in FWHM with increasing energy. 
FWHM values in the range 5-8% at 2809 keV were 
achieved. Variations in the resolution achieved can be 
attributed to changes in the contacting methods and the 
voltages applied between the different data taking peri- 
ods. It should be noted that the detectors used here do 



not have the best energy resolution possible, since for this 
first study with unknown background a very good energy 
resolution was not considered to be essential and, hence, 
cheaper detectors were used. 



V. BACKGROUND STUDIES 



To understand the observed spectrum and disentan- 
gle the individual contributions, a background model has 
been developed. All materials used were measured for 
contaminants in the LNGS Ge-detector facility, though 
some of them could only be measured after the start 
of data taking with the prototype. As a consequence 
of these measurements, the pertinax holders and lemo 
cables were replaced by delrin holders and kapton foils 
(between periods A and B). No contamination of the 
CdZnTe could be detected with the Ge facility. With 
the known activities of contaminants in the individual 
components, extensive Monte Carlo modelling based on 
GEANT4 was performed to describe the observed spec- 
trum (Fig. [3]). 

By far the largest background in the 2-3 MeV region 
evolves from the passivation paint on the detector sur- 
face. The precise prediction of this contribution varies 
slightly due to the unknown paint mass and the inhomo- 
geneous paint thickness, which affects the alpha-particle 
simulation in particular. However, there is a slight ad- 
vantage associated with the contaminated paint as the 
detector effectively acts as a self-calibrating device. The 
installation of the VME system at the start of period C 
significantly increased the timing resolution achieved for 
event read-out (from wlms to <10/is) permitting the 
observation of /3 — a coincidence events from ^^''Bi. This 
isotope originates from the Th-decay chain, present in 
the passivation paint, and beta-decays with an endpoint 
of 3.3 MeV; the daughter isotope, ^^^Po, alpha-decays 
with a half-life of 164. 3^s, releasing a 7.7MeV a. The 
rate of event pairs observed in the period C data-set was 
consistent with the measured activity of a paint sample. 
Furthermore, it shows that any possible dead layer at 
the detector surface is insignificant, otherwise the alpha- 
particles would not be detected. In the meantime, an 
alternative solution for surface passivation of the CZT 
detectors has been found and is currently being explored 
at LNGS. Initial measurements show a reduction of this 
background by at least a factor of 8, if not more, in the 
region of interest around 2.8 MeV. 

The measurements of pertinax contaminants prompted 
the replacement of all pertinax components with delrin, 
leading to a reduction of events by about a factor of 5 in 
the range 500-2000 keV, though some of this reduction 
can be attributed to the exchange of the lemo cables. 
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FIG. 3: (Color online) Summed data spectrum for period 
A, compared to simulated activities of the individual compo- 
nents. This clearly shows that events from the passivation 
paint on the detector surfaces contribute a significant back- 
ground to the range of interest (2-3 MeV). The excess of 
events below 3 MeV is thought to be due to backgrounds as- 
sociated with the lemo cables and solder. 



VI. ACHIEVABLE ENERGY RESOLUTION 

The underground studies are not currently limited by 
energy resolution and, therefore, it was not considered 
necessary to use the highest quality of crystals in this 
set-up. However, as background levels are reduced the 
energy resolution will become important, since a sharp 
peak is especially important in reducing the contribution 
of the irreducible background of 2i'Pf3 events to the 0^(3(3 
peak region. The fraction of 2i//3/3 events in the peak 
region, as a function of energy resolution (FWHM) can 
be approximated by [l0[ 



rue \ Q 



(1) 



With this in mind, additional studies were performed 
outside the underground experiment to determine the 
resolution achievable with CZT coplanar-grid detectors 
and to investigate possible improvements. 

Fig.|3]shows a ^^*Th spectrum measured with a typical 
'medium quality' CZT detector, resulting in a resolution 
of 2.6% at 2614 keV but resolutions as good as 2.1% have 
been measured with COBRA crystals. With such a res- 
olution, for the case of 90% ^^^Cd isotopic enrichment, 
2i'P(3 decays will only contribute 2 x 10"'' counts/(kg 
yr) to the 59 keV wide signal region (calculated using the 
observed half-life of 2.7 x lO^^yrs ^T|). This is already 
well below the required background levels shown in Fig. [1] 
but detectors with still better energy resolution are com- 
mercially available and are being considered for use in a 
future stage of the experiment. 

A further experimental option to improve the energy res- 
olution is cooling of the detector. This might be espe- 
cially important for searches in the low energy range, 
for decays such as two neutrino double electron capture 
(2j/ECEC) that produce a signal below 100 keV. First 
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Energy (keV) 

FIG. 4: Energy spectrum of a -^^Srpj^ calibration source ob- 
tained with a 1 cm"^ CZT detector at room temperature of 
around 24 degrees. The energy resolution (FWHM) at the 
2614 keVline is 2.6%. 



measurements of cooling from 24°C to 10°C revealed an 
improvement in energy resolution of a factor two below 
100 keV and an improvement of 5 % on the typical reso- 
lution at 2809 keV. 



VII. DATA ANALYSIS 

The data analysis consists of two independent parts: 
simulation of the possible double beta decays to deter- 
mine detection efficiencies and a maximum likelihood 
peak search. 

The predicted signals in the crystals were determined 
through a GEANT4 based Monte Carlo simulation util- 
ising calculations from the Fortran DecayO code p^ . 
Oi/PP decays to ground state (g.s.) and excited states^ 
were simulated for each of the candidate isotopes con- 
tained in natural CdZnTe. For (3~ (3~ transitions, calcu- 
lations based on the light Majorana neutrino exchange 
mechanism were used for (0^ — > 0^) transitions, whilst 
right-handed currents were used in the calculation for 
(0"^ 2+) transitions. As there is no general connection 
between ground state and excited state matrix elements, 
they must be explored separately for each isotope. 
The energy, £^pcak, and intensity of the dominant peak for 
each (3~ P~ decay were determined from these simulations 
and are given in table [ill The efficiency for observation 
of the full peak energy, e, determined from the peak in- 
tensity generally decreases with increasing peak energy. 
For decays to excited states, gamma escape probabilities 
also play a part. 
For the isotopes that decay through /3+/3+ transitions, 



^ Only excited states already included in the DecayO code were 
used, therefore some transitions, namely those to the higher ex- 
cited states of ^■^''Te, have been omitted from this analysis. 
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Isotope Decay 


-^'pcak 

(MeV) 


(%) 


e 

(%) 


Fit range 
(MeV) 


"^Cd to g.s 
i30Te to g.s 


2.809 
2.529 


4.7- 7.6 

4.8- 7.8 


66.5 
70.9 


2.21-3.20 
2.21-3.20 


"°Te to 536 keV 


1.993 


5.0-8.5 


61.2 


1.70-2.28 


"''Cd to 1294 keV 


1.511 


5.4-9.4 


74.4 


1.20-1.78 


"•^Cd to 1757 keV 


1.048 


6.0-11.0 


60.4 


0.90-1.30 


™Zn to g.s 
"*Te to g.s 


1.001 
0.868 


6.1-11.3 
6.5-12.2 


93.3 
94.8 


0.60-1.30 
0.60-1.30 


"^Cd to 2027 keV 


0.778 


6.8-12.9 


67.4 


0.50-1.20 


"''Cd to 2112 keV 


0.693 


7.1-13.8 


77.3 


0.50-1.00 


"'^Cd to 2225 keV 


0.580 


7.7-15.4 


76.6 


0.50-1.00 



TABLE II: Specifications of tfie fits performed for /3~/3~ 
peaks. Values given for each decay are the peak energy, -Epoak, 
the FWHM (Ai7peak) at that peak energy (range for all 6 data 
subsets), the efficiency, e, for observing that peak determined 
from simulations and the energy range over which the fit was 
performed. Note that decays not separated by a horizontal 
line were fit together. 




Energy (keV) 



FIG. 5: (Color online) Simulated spectrum in crystal 1 for 
^^Zn /3+EC decays to the ground state for ideal resolution 
and convolved with the energy response for data subset lA. 



double electron capture (EC/EC), single electron cap- 
ture (/3+EC) and double positron (2/3+)-transitions were 
considered when energetically possible. In general the 
predicted spectra for these decays are significantly more 
complex than those for P~ f3~ transitions, without clearly 
dominating peaks. As an example, the simulated spec- 
trum for ^''Zn /3+EC decays to the ground state is shown 
in Fig. [5l both for ideal resolution and convolved with 
the energy response for data subset lA. Although there is 
a peak at the Q-valuc for this decay (1096 keV), there are 
also single and double 511 keV escape peaks, and a peak 
at SllkeV due to gammas produced from the anihilla- 
tion of positrons produced by these decays in the other 
three crystals. 

Detailed simulations and ex-situ measurements of back- 
ground contaminants were used to characterise the mea- 



sured background continuum, as described in section [V] 
An exponential of the form y — A + Be~^'' was found to 
describe the data well at higher energies. Below 500 keV 
the fits were less satisfactory due to the presence of a 
number of low-energy gamma peaks and the 4-fold for- 
bidden decay of ^^'^Cd 13] so a 500 keV threshold was 
enforced. Thus, decay modes with no significant peaks 
above 500 keV, namely those of ^"*Cd and ^^"'Cd, were 
omitted from the analysis presented in this paper. 
A study of the residuals of the background fits showed 
some evidence for additional peaks in the continuum at 
610 keV and 1120keV, relating to gamma lines from ^^^Bi 
decay. Due to the combined efl^ects of decay branching 
ratios, energy resolution and reduced efficiency for stop- 
ping higher energy gammas in a single crystal, no other 
background gamma peaks within the analysis region were 
predicted through simulations. Therefore, two gaussian 
peaks were added to the description of the background 
with fixed mean and width determined from the relevant 
resolution function. The amplitudes of these peaks were 
treated as additional fit parameters. As a cross-check, 
fits for the exponential background were performed with- 
out these additional background peaks and in all cases 
resulted in a poorer, or negligibly different, chisquared 
probability. 

A maximum likelihood fit was performed to determine 
the most likely number of signal events, 6s, over the com- 
bined data-set. Parameters describing the background 
were allowed to vary between crystals and data collec- 
tion periods, but the Oz^/?/3 signal rate was assumed to be 
constant throughout, i.e. Different background param- 
eters were applied to different data subsets to allow for 
the varying background rates (indicated in tabic |T| but 
the normalised background distributions fitted to each 
data subset were found to agree within errors for each fit 
scenario. 

For (3~ l3~ modes, ds enters the fit through the ampli- 
tude of a gaussian peak with width determined by the 
calibrated resolution of the relevant data subset. The 
range of peak widths (FWHM) for each fitted peak are 
given in table[TTlalong with the energy range used for each 
fit. Simulations showed that a range of (iJpoak ± 3AE) 
or greater was required for each peak-search in order to 
adequately characterise the background continuum. The 
close proximity of some predicted signal peaks required 
these ground-state transition signals to be determined si- 
multaneously: ^^^Cd and -'^^^Te were fit together and the 
^•^Zn and ^^^Te peaks were also fit simultaneously. Tran- 
sitions to excited states are expected to be significantly 
suppressed with respect to ground state transitions due 
to phase-space arguments, so any contribution of excited 
state decays to the fitted peaks for ground-state transi- 
tions is assumed to be negligible. The limit for each signal 
arising from a transition to an excited state was deter- 
mined in a separate fit. The high-threshold (> 600 keV) 
data set was used for all /3~/3~-mode peak searches ex- 
cept the decays to the third and fourth excited states of 

"6Cd. 
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The majority of spectra predicted for /?+/3+-mode decays 
have multiple peaks, each significantly smaller in ampli- 
tude than those arising from (3~P~ decays to the ground 
state, thus justifying separate treatment in the analysis. 
Due to the complexity of these spectra, to determine the 
most likely number of signal events, 9s, the most likely 
scaling factor for the entire simulated spectrum, was ex- 
tracted from the likelihood fit. For each decay the sim- 
ulated spectrum was normalised to unity and convolved 
with the relevant resolution function for each data sub- 
set. The range for each fit (as given in table IIV|1 was 
selected to include all the dominant peaks. For tran- 
sitions where the simulated spectra showed a dominant 
peak in the region 500-600 keV, the low-threshold data 
set was used; in all other cases, the high-threshold data 
set was used. 

For all P~ (3~ and modes analysed, a 90% confi- 

dence limit on the half-life, Thaif, was determined from 
the fitted number of signal events, 0s ± (5^, where 6s is 
one sigma uncertainty in the fit, for each decay under 
investigation. 

_ iViso^iivceln2 

" + 1.285,) ^ ' 

Here A^igo is the number of candidate nuclei per crystal 
for the given decay, tuvc is the total duration of data 
collection in crystal-years, and e is an efficiency factor 
determined from simulations. For /3~/3~-searches, e is 
the fraction of simulated events in the peak region (see 
table HI]) . However, for /3+/3+-searches e = 1 since the 
whole simulated spectrum, normalised to unity, is used 
in the fit. 

For each fit, a chisquared goodness of fit test was per- 
formed. However, due to the low statistics, this param- 
eter was not expected to follow a true chisquared dis- 
tribution. Therefore, the distribution of the statistic 
was determined by Monte Carlo for each fit in order to 
calculate the fit probability. The x^i a-nd its respective 
probability, determined for each /3~/3~-fit is included in 
table IHll For the /3+/3+-fits the probabilities for each fit 
were also all >80%. As a cross-check, fits were repeated 
with 6s fixed to zero for each signal, all of which resulted 
in either negligible change or a decrease in the goodness 
of fit. 

A detailed study of possible systematic effects was per- 
formed and the dominant uncertainties were found to be 
those that affect the number of candidate nuclei. Un- 
certainties in energy resolution and livetime, and possi- 
ble biases in the fit procedure, were all found to have a 
negligible effect on the analysis. However, the possible 
existence of a dead-layer at the surface of the crystals 
could reduce the active volume by up to 10%. Observa- 
tions of a-particles from the passivation paint indicate 
that the dead layer is probably smaller than this but the 
effect was taken into account in a conservative manner 
by using iViso x 0.9 in the limit calculation. Due to the 
production process, the zinc content is only known to be 
in the range 7-11% resulting in an uncertainty in both 
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FIG. 6: (color online) The fitted 90% confidence limits for 
"''Cd and ^^"Te Ov(3f3 decay events (top) and ^^°Te Oi/EC/EC 
to ground state events (bottom) with the total high-threshold 
data set shown for the range of each fit. Also shown separately 
is the contribution from the fitted background. 

the number of zinc nuclei and the number of cadmium 
nuclei. To ensure conservative half-life limits, 7% zinc 
content was used when calculating the number of zinc 
nuclei, and 11% zinc content was used in calculations for 
cadmium isotopes. 



VIII. RESULTS 

TableHnishows all the p-p-decay half-hfe hmits (90% 
C.L.) calculated in this work and table lIVl shows the lim- 
its calculated for /3+/3+ decays. The limits obtained from 
the combined fit to the "^Cd and ^^^Te Oi^PP decays to 
ground state and for the ^^°Te double electron capture 
decay to ground state are shown in Fig. [SI Due to the 
small detector mass, at present the searches for ^^^Cd 
and ^■^"Te cannot compete with other running large scale 
experiments. However, half-life limits obtained for ^"^Zn 
and ^^"Te improve on existing measurements. 



IX. SUMMARY 

A new double beta decay experiment, COBRA, is 
planned using a large amount of CZT semiconductor de- 
tectors. A low rate of background events in the peak 
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Isotope and Decay 


xVdof 


P 


Ti/2 li 
This work 


mit (years) 
World Best 


"'^Cd to g.s 
^=*°Te to g.s 


1236/2950 
1236/2950 


0.99 
0.99 


3.14x10" 
9.92x10" 


1.7x10^^ [14] 
LSxlO^"* [15] 


^''"Te to 536 keV 


726/1703 


0.99 


3.73x10" 


9.7x10^2 [16] 


"^Cd to 1294 keV 


676/1703 


1.00 


4.92x10^** 


2.9x10^^ [14] 


"«Cd to 1757 keV 


615/1479 


0.99 


9.13x10^® 


1.4x10^2 [M] 


''°Zn to g.s 
^^®Te to g.s 


859/2078 
859/2078 


1.00 
1.00 


2.24x10^^ 
5.38x10" 


9.0x10^'^ [17] 
1.1x10^^ [18] 


"^Cd to 2027 keV 


756/1732 


0.96 


1.37x10" 


2.1x10^^ \M 


"'^Cd to 2112 keV 


545/1233 


0.87 


1.08x10" 


6.0x10^^ [14] 


"^Cd to 2225 keV 


545/1233 


0.87 


9.46x10^* 


LOxlO^""!" [20| 



region and good energy resolution are crucial aspects of 
the design of such an experiment. To develop this new 
approach, CZT semiconductor detectors have been op- 
erated deep underground, for the first time, to study 
their background. Using a small prototype, a background 
model has been developed and a major background com- 
ponent in the form of a passivation paint on the detector 
surfaces was identified. Alternatives are now available 
and major improvements are expected soon. 
Studies of the attainable energy resolution showed that 
the contribution of the irreducible background of 2i'f}(3 
can be kept to a negligible level. A 4.34kg-day data set, 
collected with four 1 cm'^ crystals, has been analysed to 
determine limits on various neutrinoless double beta de- 
cay modes of seven different isotopes. Despite the small 
detector mass, only 26 g of CZT, these data have yielded 
four improved half-life limits for decays of ^^Zn and ^^"Te. 
After submission of this, paper new limits for ^^°Te have 
been presented [23 |. 

In the near future, 64 CZT detectors will be running and 
improvements on all the limits presented in this paper 
can be expected. In addition to the increased detector 
mass, new criteria based on coincident energy deposits in 
time and space will allow better rejection of backgrounds, 
and a full characterisation of each individual crystal will 
help reduce systematic uncertainties. 
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TABLE III: Results from fits for decay peaks. The 

X^/DoF goodness of fit parameter, and respective probabil- 
ity determined through Monte Carlo are included. The 90% 
confidence limits have conservative systematic uncertainties 
applied and are compared to the world best limits. ^ Quoted 
limit is 68% not 90%. 



Isotope and Decay 


Range 


Ti/2 limit (years) 




(MeV) 


This work 


World Best 


^^Zn 0uf3+EC to g.s. 


0.5-1.3 


2.78x10^'^ 


2.4x10^** [17] 


'^'^Zn 0i/2EC to g.s. 


0.7-1.3 


1.19x10^^ 


7.0x10^'' 117} 


^^°Te Oi^/3+EC to g.s. 


0.5-2.0 


1.21x10^^ 


2.2x10^** [2JJ 


^2°Te 01/2EC to g.s. 


1.2-2.0 


2.68x10^''' 




^2°Te 0i^2EC to 1171keV 


0.5-2.0 


9.72x10^^ 




i^^Cd 0uf3+f3+ to g.s. 


0.5-2.0 


4.50x10" 


2.4x10^" [22] 


i"^Cd 0uf3+EC to g.s. 


1.4-3.0 


7.31x10^** 


3.7x10^" [22] 


^°^Cd 0i^2EC to g.s. 


1.4-3.0 


5.70x10^^^ 


1.5x10" [23] 


i^^Cd 0uf3+f3+ to 512keV 


0.5-2.0 


1.81x10" 


1.6x10^° [22] 


^'''^Cd 0^/3+ EC to 512keV 


0.8-2.0 


9.86x10" 


2.6x10^" [22] 



TABLE IV: 90% confidence limits obtained for all de- 
cays analysed in this work with conservative systematic un- 
certainties applied, compared to the world best limits. New 
world best values from this work are shown in bold. The 
energy range used for each fit is also included. 
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